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Abstract 

The issues of scaling symmetry and critical point behavior are studied for fluctua- 
tions about extremal charged black holes. We consider the scattering and capture of the 
spherically symmetric mode of a charged, massive test field on the background spacetime 
of a black hole with charge Q and mass M. The spacetime geometry near the horizon of 
a \Q\ = M black hole has a scaling symmetry, which is absent if \Q\ < M, a scale being 
introduced by the surface gravity. We show that this symmetry leads to the existence of 
a self-similiar solution for the charged field near the horizon, and further, that there is a 
one parameter family of discretely self-similiar solutions . The scaling symmetry, or lack 
thereof, also shows up in correlation length scales, defined in terms of the rate at which 
the influence of an external source coupled to the field dies off. It is shown by constructing 
the Greens functions, that an external source has a long range influence on the extremal 
background, compared to a correlation length scale which falls off exponentially fast in 
the \Q\ < M case. Finally it is shown that in the limit of A = (1 — ^p-) 5 — > in the 
background spacetime, that infinitesimal changes in the black hole area vary like As. 
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1. Introduction 



Recently, numerical studies of gravitational collapse, have shown scaling and critical 
point type behavior in the formation of zero mass, neutral blackholes [0J , [p] . These studies 
showed several interesting properties; in the zero mass limit, the wave form of the collapsing 
wave always evolved to a particular form, which was discretely self-similiar in appropriate 
time and space variables. The mass of the formed black hole had non-analytic dependence 
on a variety of parameters measuring the difference in the strength of the wave from some 
critical value, and the exponent was found to be universal. 

Suppose that charged particles collapse to form a black hole. In this case, the black 
hole must have a mass which is greater than or equal to its charge. Would critical point 
type behavior be seen in fluctuations about the minimal area? Or more generally, would 
such behavior be seen in the interactions between charged wave packets and an already 
existing charged black hole, in the extremal limit? 

There are two geometrical reasons why this might occur. First, the spacetime geome- 
try near the horizon of a Q = M black hole has an infinite throat (the metric approaches 
a Robinson-Bertotti metric). The throat has no scale, and the metric has a dilatation 
symmetry, which means that test fields on this background will have a scale invariance. 
By contrast, the geometry near the horizon of a non-extremal black hole has a scale set by 
the surface gravity k. Second, dust with mass density equal to its charge density can be 
placed in arbitary configurations and will stay in equilibrium with other configurations of 
such dust , and with arbitrary distributions of charge equal to mass black holes. There is 
no particular size of charge equal to mass dust that is needed for a force balance. This is 
reminiscent of the picture of fluctuations on all length scales occuring at a critical point. 

In this paper we will focus on the questions of scaling invariance and self-similiar 
solutions, correlation length scales, and how fluctuations in the area of the black hole 



depend on A = Jl- fa, as A -> 0. Of course one would like to have exact solutions 
describing wave packets of charged fields scattering off a charged black hole, analogous to 
the numerical work 0,0. Here, in order to make some progress analytically, we will study 
a charged, massive, test field scattering off a fixed black hole background with charge Q and 
mass M. This is a consistent approach since (1) we are interested in the limit where the 
change in the area is infinitesimal, and (2) because there is already a black hole present to 
do perturbation theory around, unlike the neutral black hole case. We imagine an initial 
wavepacket which heads towards the black hole, part of which is scattered and part of 
which is captured. One is interested in the form of the captured wave, in particular to 
see if it shows scaling behavior when the background spacetime approaches extremality. 
We will first show that near the horizon of an extremal black hole, the wave does have a 
scaling symmetry when the background spacetime is extremal. Amongst these solutions 
there is one which is self-similiar, and has a translation invariance in logarithmic time and 
logarithmic radial coordinates. Further, we will show that near the horizon, there exists a 
set of eigenfunctions of the wave equation which each of which has a discrete self-similiarity. 

Second, addressing the issue of correlation lengths is a bit confusing- this is not (at 
least apparently) a statistical system with degrees of freedom to average over. However, 
one can construct the Greens functions for the wave equations, and these tell what the 
response of the test field is to a source. The Greens functions show that on the extremal 
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background the influence of a source is long range, whereas on a background with Q < M, 
it falls off exponentially fast. Thirdly, the captured part of the wave adds mass and charge 
to the black hole. We will show that the resulting change in the area of the black hole goes 
like A 2 in the limit where A — > 0. 

To avoid repeated absolute value signs, we will take the charge of the black hole to be 
positive. We will use units with G = 1. 

2. Description of the System 

Consider a charged scalar field on a Reissner-Nordstrom background spacetime, with 
metric given by 

ds 2 = X 2 (-dt 2 + dy 2 ) + R 2 dQ 2 , X 2 = 1 - ^ + % (1) 

R R 

Here y is the usual tortoise coordinate, dy = dR/X 2 , R — R(y) and Q an M are the electric 
charge and mass of the spacetime. Assume that the charge Q is positive. Choose the gauge 
of the electrostatic potential such that it vanishes on the horizon, R = Rh, 

At = o -§, O = ^- (2) 

it jxh 

The quantity (fi is the then the difference in the electrostatic potential between the horizon 
and infinity. In the context of the first law of black hole thermodynamics, (j) is conjugate 
to the electric charge Q. Note that for Q = M, 

M 

A t = X = l-—, <f> = l. (3) 

The matter action, for a scalar field of mass m and charge e, is given by 

S m = J ^(D a $D a <f>* - m 2 $$*), (4) 

where the gauge covariant derivative -D $ = (d a — ieA a )Q. The equation of motion for 
the scalar field is 

V a V a $ - 2zeA a V a $ - (m 2 + e 2 A a A a )§ = (5) 

We will consider spherically symmetric waves on the background ([!]). Let ip = 
Then the equation of motion for the field (|5|) becomes 

d 2 y il> - d 2 tP + 2ieA t d t i> + (e 2 A 2 - m 2 X 2 - V grav )ijj = 0, (6) 

with the potential V grav given by 

_d 2 R_2M Q 2 

Vgrav ~ "IT ~ ~W X (1 ~ MR ] (7) 
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Lastly, we summarize the behavior of the system near the horizon. In the tortoise 
coordinates, the black hole horizon, R = Rh, is at y — > — oo. One finds near the horizon 
that: 

for Q = M, 

M 2 M 
R-M -> , and A t = A , (8) 

V V 

while for Q < M, 

R — Rh — > R H e 2Ky , A 2 -> 2nR H e 2Ky , A t -> c/> e 2Ky , (9) 

where /c = ^9fjA 2 |^ H is the surface gravity at the horizon , which vanishes for the extremal 
black hole. 

3. Scaling Behavior and Self-Similiar Solutions near the Horizon of Extremal 
Black Holes 

From the asymptotic behaviors of the gauge potential and metric function, it follows 
that the wave equation (|6|) for the scalar field is invariant under the rescaling 

y — > ay, t — > at, a = constant. (10) 

This implies that there are solutions to (|J) which are functions only of the ratio t/y. 
Equivalently, in term of logarithmic coordinates i = ln£, y = ln(— y), there are solutions 
of the form ijj = F(t — y). These solutions have the translation invariance, corresponding 
to self-similiarity in the original t,y coordinates, ip(i + D,y + D) = tf;(i,y) for any D. 
By contrast, the field equation on the nonextremal black hole spacetime does not have a 
scaling invariance-the surface gravity k introduces a scale. 

The scaling invariance is a reflection of an additional dilatation symmetry of the 
metric near the horizon, ds 2 — > ^-(~dt 2 + dy 2 ) + M 2 dQ 2 , with the dilatation killing 
vector £ a = t{Mr) a + y(jr) a ■ The symmetry is actually best seen in a slightly different set 



>, in which tl 
the wave equation (W) becomes 



of coordinates, in which the metric function is used as a coordinate. Let x = 1 — Then 



-d t 2 $ + 2iexd t $ + x 2 K d x {x 2 d x <$>) + x 2 (e 2 - m 2 )$ - 0. (11) 

As x —>■ the horizon is approached and the differential equation becomes invariant under 
transformations of the form x -x,t — > at, with a = constant. The utility of the 
coordinate x is that it is consistent to include the term (e 2 — m 2 )$ . In the tortoise 
coordinates above it is not clear that one can legitimately retain this term, while ignoring 
higher order terms in the inversion relation between y and R. 
Now let w = xt and look for solutions of the form 

${x,t)=t iv F v (w) (12) 
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In terms of the logarithmic time i introduced above, the prefactor is e lut , and so the 
solution will have a discrete self-similiarity t — > t + — and x = —Inx — > x + — For the 
eigenvalue v = this is the continuously self-similiar solution discussed above, depending 
only on w. 

The wave equation (|TT|) becomes, for x <C 1, 



, w , w// , , u> iu s , 9 9 2ez^ , , 

(! - Ml)^ H + 2 H e - M2 + + (™ - e + — - -^)^ = ( 13 ) 

Analysing the solutions to this equation is a topic for future study, however, we have 
shown that the wave equation for a massive charged field on the background spacetime of 
an extremal black hole, has discretely self-similiar solutions, as the horizon is approached. 
It is worth noting that for v = the differential equation becomes 



((1 " Jfi)f» ~ 2ie f) - ( e2 - m2 )f = 0- (14) 



d f w 
dw 



For e = m it is simple to find the solution, 



f( w ) = C 1+ C 2 (§±^)- M . (15) 
M — w 

In terms of the tortoise coordinate y this is 

$ = Cl + ° 2{ T~^t/tj rM > as y (16) 

A general solution to the wave equation near the horizon can be written as a sum of the 
eigenfunctions (|T2|). However, an arbitrary sum will no longer be discretely self-similiar. 
So, perhaps the most interesting question is, do generic wave packets starting in the flat 
region (or ANY packets for that matter), evolve into a packet which is a special sum of 
the modes ([12]), such that the sum is discretely (or continuously ) self-similiar? We do 
not yet know the answer to this, but one can imagine at least two ways in which this could 
happen. When the eigenvalue problem is solved with appropriate boundary conditions for 
a wave packet incoming from the flat region, it may be that the eigenvalues v are actually 
quantized , say v n = rijr. Then the solution would have a discrete self-similiarity with 
D = . Alternatively, it could be that the imaginary part of the frequency v is positive, 
so that the lowest frequency dominates at late times. Or, it may be that the evolution 
from wave packets "at infinity" is not self-similiar! It will certainly be of interest to resolve 
this question. 

4. Scattering and Asymptotic Solutions 

The scaling symmetry can be displayed in terms of the greens function for the scalar 
field equation. The greens function describes how the field propagates in response to an 
external source. We will show that the greens function for the Q = M case has long 
range 1/y correlations, compared to the greens function for the Q < M case, in which 
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correlations die off exponentially with scale k~ x . To this end, we will find eigenmodes of 
the wave equation in the asymptotic regimes, and use the solutions to construct the greens 
functions. Further, the scattering behavior of the eigenmodes will be anlyzed, to determine 
what scatters and what is captured. From this, the change in the horizon area 8 A will be 
found, when a small amount of mass and charge is captured. 5 A has non-analytic behavior 
as extremality is approached, whereas away from Q = M, 5 A is linear in the added mass 
and charge. 

To put the wave equation in the form of a scattering problem, first, fourier transform 
in time; let 

dtl p = _ e( p o )^ (17) 

Then @ becomes 

dfo + [(u - e<f> ) 2 + 2e(uj - e<j) )A t + e 2 A 2 t - m 2 X 2 - V grav ]^ = (18) 

or, 

d 2 i; +[k 2 -V- V grav ]ijj = , k 2 + m 2 = u 2 (19) 

where 

V -> \{ewQ - m 2 M) - (e 2 - m 2 )-^, R -> oo (20) 

At large distances R —>■ y and the form of the potential V is the same in all cases. V grav falls 
off like R~ 3 . However, near the horizon there is a qualitative difference for the extremal 
and non-extremal cases: 
For Q < M, 

V -> 2ue(j> - e 2 (p 2 - m 2 - (2e<p (uj - e<j> ) - m 2 2KR H )e 2Ky + 0(e 4Ky ) , y -> -oo. (21) 
For Q = M, 

M M 2 

V -> -m 2 - e 2 + 2euj + 2e(u - e) (e 2 - m 2 )^, y -> -oo. (22) 

y y 

Analysis of the potentials shows the following qualitative features of the scattering 
problem: At large y, the potential falls off quite slowly (like y~ 1 )- A WKB approximation, 
shows that the transmission is exponentially suppressed if the incident wave is under 
the barrier. For the purposes of the following discussion then, it will be sufficient to 
approximate the capture cross section as a step-if the wave is over the barrier, it is captured, 
and if it is under the barrier, the wave is scattered. Essentially, we are working in a 
geometries optics approximation. (However, we know that the approximation is a good 
one here, because the scattering problem is similiar to that in || , [|4j , in which the scattering 
is worked out in analytic and numerical detail.) 

So what we need is the criterion for an eigenmode to be over the barrier. Now, this is 
not quite a standard scattering equation, because the height of the potential depends on 
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the incident wave frequency u. But studying the potential, one finds that to be over the 
barrier, a wave with frequency 10 must satisfy 



■m 2 A -I- f 

(u-e<p o )> -=Z±= ,e>0 (23) 

e vl - A 2 

where A = ^1 - (24) 

A is a parameter which measures how close the spacetime is to extremal, e is any number 
greater than zero, and merely insures that a wave packet centered on the frequency ui 
reaches the horizon in finite time. This condition is the same as one finds from analyzing 
paths of charged particles, which is a much simpler way to see the results. 

The scaling symmetry discussed above shows up in the form of the solutions near the 
horizon, and also in correlation lengths. Of course, this is not a statistical or quantum 
mechanical system (though the fact that Hawking radiation makes it seem like one is 
intriguing) so to find correlation lengths we can't take an ensemble average. However, we 
can compute the Greens function for the wave equation, which tells how the classical field 
evolves in response to a general source. Quantum mechanically, the two point correlation 
function is the greens function. Next we show via the Green's function, that the influence 
of a source dies off expontially for Q < M, whereas there is a long range tail for Q = M. 

Let a wave be incident on the black hole, with ip ~ _i e -»(w-e<£ )t—sfcfl ag _^ ^ 



with (u - e(p ) satisfying (||). Then as y ->• -oo, i/j -> \J ( UJ - e $ ) e ~ i{ "~ e4 ' o){t+v) ■ The 
normalization follows because the wronskian of ([TSp is constant, and using the fact that 
the amplitude of the captured wave is much greater than the amplitude of the scattered 
part. This is only the leading term in the asymptotic solution for ip. More information will 
be needed, so next we find the leading nontrivial behavior of the wave near the horizon. 
For Q = M and for y <C — M, the wave equation becomes 



, 2 . M 
(u> — e) — 2e(uj — e) — 

y . 



tp = 0. (25) 



The inward propagating solution correct through order ^ is 

= e ~ lb e 2 ^-^y (26) 



e 2 M 2 1 

where S = (u - e)y + eMln(-y) + — --. (27) 

2{u-e)y 

The outward propagating mode is given by 

• o _ eM 

ip 2 = e lb e »(«-—)» (28) 
For Q < M and ny —1, the wave equation becomes 

dfa + [(u - e<p o ) 2 + P^] ip = (29) 
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where = 2e<p (u! — e<j) ) — 2kw?Rh. must be positive for the wave to actually get 
over the barrier, by (|23|) . There are two (potentially) small parameters here, k which goes 
to zero in the extremal limit, and {u — e<p ) which we want small to be adding a small 
amount of mass to the black hole. The regime of interest will be k <C {uj — ecf) ), to approach 
the extremal black hole, and this is included in the P w > range. 

The inward propagating solution, correct through terms of order e 2liy is 

where A = T , f Y2 , 2 ( 2 ^" e ^) + i). For k < (u - e<j> ), A -> ^ - ™ 2r " « + 

^ 2{w-l^> ) — ^ ou ^ war< ^ propagating mode is the complex conjugate, 

^4=^3- (31) 

As a consistency check, the wronskian of each of the above solutions ipi(y) is constant, as 
it should be. 

The (advanced) Greens function can now be constructed. For Q = M this is 



G(t, y; t', y') = - -e(At)S(-Ay)S(At + Ay)e- ieMlog (v/v') 
2 



-hi- ieM )( y --i) 
2 y 



+ e -^(A*+Ay) e I(l-ieM)(^-l) 



- -e(At)Q(Ay)Q(At - Ay)e ieMl ° 9i V^ 



* y 



For example, the field response to a delta function source at y Q , t is 



(32) 



1>(y,t) « e - ieM ^(^)[i(l - zeM)(^ - 1) - i e -i(i-«-^ e -*i£(**+*v)(i + (33) 

4 y 2 2 y 

The solution shows the scaling symmetry, and long-range correlations, i.e., the effect of 
the source falls off like y~ l . G (or if)) has a piece which is independent of time, and a piece 
which goes to free oscillations at a frequency eM/y a , which depends on the location of the 
source point. 

By contrast, for Q < M there is no scaling symmetry, and the influence of the source 
falls off exponentially fast, like e 2Ky . The green's function has pieces which are oscillations 
at two frequencies, a = e4> e AKy , and kj, where 7 = m f~ H The first depends on the 
location of the source point, and the second goes to zero in the extremal limit. Let 
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h(Ay) = e 2hiAy — 1. Then the greens function for Q < M is 
G(y, t; y', t') =Q(At)Q(-Ay)Q(At + Ay^e"^^ 



^ ie i f i(At+Ay) e ±(l-i~/)h(Ay) 



^(At-Ay) e $(l+i-y)h(Ay) 



• (34) 



_^l e -i^(At+Ay) e -t(^±llh(Ay) 

y 2 

+ e(At)6(Ay)e(At - Ay)\e l ^ h{Ay) 

_^l e -i^(At+Ay) e -i^lh(Ay) 

For example, the field configuration at large negative values of y, due to a point source, can 
be read off from the first two lines, showing that ift approaches free oscillations exponentially 
fast. 

5. Critical Point Behavior? 

For formation of a zero mass, neutral black hole, the numerical studies looked at how 

the mass of the new black hole varied with the parameters of the incident wave. It was 
found that this behavior was non-analytic (and universal). When forming a charged black 
hole, the relevant quantity may be fluctuations about the minimal area. In the present 
context, let us look at infinitesimal changes in the area of an already existing black hole, in 
the limit where the black hole approaches extremal. (This is the analogue of the approach 
to the putative critical point.) For a black hole with general charge and mass, the horizon 
radius is Rh = M(l + A) and the area is AnR 2 H . If small amount of mass 5M and charge 
5Q are added, the change in the horizon radius is 



SR H = SM + M ] /A 2 + 2 5 -^-2^- MA, (35) 

where terms of order 5M 2 ,5Q 2 have been neglected. For a wave carrying 5Q = e to be 
captured, as discussed in (|2^) , it must have 

2 

777 

SM = u = ecf> + — (A + e), (36) 

e 

where e — > to add the minimal possible mass. Now there are two cases. If one fixes A 
and then considers 5M, 5Q A, then as expected, one finds a formulae for the change in 
radius which is linear in the perturbations to the mass and charge, 

6R H « e(</> -^ K )+ m ^(A + e)(l + I). (37) 

On the other hand, for the case of interest here, m and e are still small compared to M 
and Q, but in addition, A — > 0. Precisely, for A <ti m 2 /(Me), 



A7? w «r + A /?^[(l-4)A + e]i. (38) 



rrr 
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Therefore, in the limit where the change in the area is as small as possible (e — ► 0), the 
variation in the horizon area has non-analytic behavior, as the extremal background is 
approached. Of course, ( |3"8D could be written as linear md, where 5 = (1 — ^r)*- Here 
A appeared to be a natural choice because the horizon area is polynomial in A. The 
suggestion in (|38| ) is that tuning through the background spacetimes as A — > is like 
tuning the magnetic field to a critical value. 

The First Law states that 5M = -^SA + (p SQ. k and </> play the roles of the 
temperature and an electric (or chemical) potential, so derivatives of SRh with respect to 
K, 4> are also of interest. Instead of using M and A as the two independent variables to 
describe the state of the system, we switch to the variables k = M ^ A ^ and <p = i^-- 
Then, for example, the analogue of the specific heat is 



,'dSA\ A „ fd5R H \ n „ 2Mm 2 r/i e 2 . . 

=4. K R H (—) =-2.^ — [ ( l--)A + e].. (39) 



6. Discussion 



Consider the class of spacetimes consisting of a charged black hole, parameterized by 
M and A, interacting with charged matter. We have been looking at various properties of 
this system, that suggests that the point A = is like a critical point. To examine this, 
we move away from this point (look at Q < M spacetimes ) and probe the system with 
charged test fields, to see the behavior as A — > 0. It was seen that fluctuations in the 
area of the black hole have non-analytic behavior in this limit. We showed that correlation 
lengths, defined in terms of the classical Greens functions, are long range on the A = 
background, and decay exponentially for spacetimes with A > 0. Further, we showed that 
spherically symmetric packets of the test field evolve to configurations which have a scaling 
symmetry near the horizon, if and only if the background has A = 0. Near the horizon, 
the eigenmodes can be chosen such that each mode has a discrete self-similiarity. This is 
interesting, because the numerical studies of formation of neutral black holes showed that 
the collapsing field was discretely self-similiar, near the critical point. In the present study 
we don't know if the same is true; does an incident wave packet evolve into a special sum 
of the eigenmodes, such that the sum has a discrete self-similiarity. This is an interesting 
open question. 

Fluctuations about zero mass is a limiting case of fluctuations about the minimal area 
(irreducible mass) , when the black hole is neutral. This would suggest that a key feature 
to criticality is extremality. However, there is another possibility which is interesting to 
think about. A black hole with Q = M in a spacetime with positive cosmological constant 
is not extremal. However, it does have the property that the surface gravity of the black 
hole is equal in magnitude to the surface gravity of the deSitter Cauchy horizon, i.e., the 
two temperatures are the same. Geometrically, the spacetime geometry has an infinite 
throat near the horizon of the black hole, similiar to the geometry of the throat discussed 
here. Therefore, one might expect that the behavior of a charged test field would be the 
same as in the present case. If this is true, then the key ingredient would be equality of 
the Hawking temperature with the background. 
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